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Based  on  the  assumptions  of  inviscid  and  incompressible  fluid,  irrotadonal  flow,  and 
infinitesimal  wave  amplitude,  Stokes  (1846)  found  a  solution  to  the  water-wave  problem  with  a 
uniformly  sloping  impermeable  boundary  (beach).  The  solution,  often  termed  the  Stokes-mode 
edge  wave,  can  be  written  in  terms  of  velocity  potential,  $,  as 


<j>  (x,y,z,t)  =^sinpe  ky°0*P  +  k2sinP  sin  (kx-rnt). 


0) 


where  A  is  the  amplitude  of  wave  runup  distance  along  the  beach  surface,  co  is  the  wave  angular 
frequency,  k  is  the  wave  number  in  the  longshore  direction,  f)  is  the  beach  slope  from  the 
horizontal,  and  the  coordinates  (x,  y,  z)  point  to  the  alongshore,  offshore,  and  vertically  upward 
directions,  respectively,  as  shown  in  Fig.  1.  Equation  (1)  indicates  that  edge  waves  propagate 
parallel  to  the  shoreline  with  their  crests  pointing  offshore.  The  maximum  wave  amplitude  occurs 
at  the  shoreline,  y  =  0,  and  decays  exponentially  offshore  with  an  e-folding  distance  of  (k  cos  P)1. 
The  corresponding  dispersion  relation  is  given  by 

co2  =  gk  sin  p.  (2) 

This  dispersion  relation  indicates  that  edge  waves  are  dispersive,  i.e.  the  group  velocity,  9co/3k,  is 
a  function  of  k.  Rotating  the  coordinates  from  the  horizontal  about  the  longshore  axis  by  P,  (1) 
and  (2)  can  be  expressed  by 

<j>  (x,Y,t)  =  —  e' kY  sin  (kx  -  cot),  (3) 

co 


and 
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co2  =  yk, 


(4) 

where  Y  is  the  coordinate  pointing  offshore  along  the  beach  surface  as  shown  in  Fig.  1,  and  y  =  g 
sin  P  represents  the  down-slope  component  of  gravity.  With  the  rotation  of  coordinates,  the 
Stokes-mode  edge  waves,  (3),  are  no  longer  dependent  on  the  coordinate  Z  perpendicular  to  the  x- 
Y  plane.  Note  that  (3)  and  (4)  have  exactly  the  same  forms  as  the  linear  solution  of  deep-water 
waves.  It  was  shown  that  nonlinear  behaviors  of  edge  waves  are  also  analogous  to  those  of  the 
two-dimensional  deep-water  waves  (Yeh,  1987). 

Just  like  the  solution  of  deep-water  waves,  nonlinear  correction  in  $  does  not  arise  from  the 
next  higher-order  solution  of  the  Stokes-mode  edge  waves  (Whitham,  1976).  Nevertheless  the 
nonlinear  correction  appears  in  the  vertical  surface  displacement  tj.  Substituting  the  linear  solution 
(1)  into  the  dynamic  boundary  condition  along  the  free  surface  yields 

T|  (x,y,t)  =  A  sin  p|^  cos(kx-cot)  -  Ak  e  * 2  ^  cos  P]  +  o(A3k2;  p3),  (5) 


where  the  beach  slope  P  is  assumed  to  be  small.  The  second  term  in  parenthesis  in  (5)  represents 
the  set-down  of  mean  water  level  of  which  effect  decays  exponentially  offshore  with  an  e-folding 
distance  of  (2  k  cos  P)1.  The  set-down  effect  of  edge  waves  was  verified  in  the  laboratory 
experiments  by  Yeh  (1986).  While  the  vertical  water  surface  displacement  has  the  characteristics 
represented  by  (5),  it  is  shown  in  this  paper  that  the  shoreline  profile  along  the  beach  surface  does 
not  involve  the  set-down  but  instead  the  second-harmonic  component  appears  at  the  same  order  of 
the  solution. 

Instead  of  the  standard  Eulerian  formulation,  edge-wave  motions  are  first  formulated  in  the 
Lagrangian  coordinates  with  one  of  the  axes  pointing  along  the  beach  surface  as  shown  in  Fig.  1. 
Assuming  an  incompressible  fluid,  the  equation  of  mass  conservation  can  be  written  as 


B(x,  Y,  Z)  = 
a  (a,  b,  c) 

u  u  a< 

where  the  operator 

d( 

fluid  particle:  the  coordinates  (a,  b,  c)  are  often  taken  to  be  the  initial  Cartesian  coordinates  of  the 
fluid  particle  at  t  =  0,  although  such  an  assignment  is  not  necessary.  Assuming  the  inviscid  fluid, 
the  equations  of  motion  in  the  (x,  Y,  Z)  directions  are,  respectively. 


(6) 


is  the  Jacobian,  and  (a,  b,  c)  are  the  Lagrangian  coordinates  to  identify  a 
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(7) 


dh 

1  3(P.Y.Z) 

at2 

p  3  (a,  b,  c) 

a2v 

i  a  (x,  p,  z) 

at2  “ 

p  a  (a,  b,  c) 

a2z 

i  a(x,  y,  p) 

at2 

p  a  (a,  b,  c) 

g  sin  (3 , 
g  cos  p , 


with  the  boundary  conditions. 


Z  =  0  at  c  =  0, 
p  =  0  at  c  =  b  tan  p, 


(8) 


where  p  is  the  pressure.  An  advantage  of  the  Lagrangian  formulations,  (6)  and  (7),  is  that  a 
moving  runup  waterline  can  be  described  with  the  independent  variables  (a,  b,  c),  hence  the 
location  of  the  waterline  is  a  known  quantity  a  priori. 

Using  the  straightforward  (regular)  perturbation  method  with  the  small  parameter  e  = 
0(tan2P)  =  0(Ak),  the  second-order  solution  of  the  Stokes  mode  in  the  Lagrangian  coordinate 
system  is  found  to  be 

x  =  a  -  A  e  ‘ kb  sin  (ka  -  cot)  +  0(e3), 

Y  =  b  -  A  e  - kb  cos  (ka  -  cot)  - 1  A2k  e  * 2  kb  +  0(e3),  (9, 

Z  =  c  +  0(£3), 


where  A  is  the  wave  runup  amplitude  along  the  beach  surface.  The  shoreline  profile  is  found  by 
evaluating  (9)  at  b  =  0  and  c  =  0: 


a  -  A  sin  (ka  -  cot)  +  0(e3), 

(10) 

-  A  cos  (ka  -  cot)  -  j  A2k  +  0(£3). 

(11) 

The  shoreline  profile  in  the  Eulerian-coordinate  system  can  be  retrieved  by  substituting  (10)  of  the 
form. 


a  =  x  +  A  sin  (ka  -  cot)  +  0(e3), 


(12) 


3 


into  (1 1).  Using  some  trigonometric  identities  and  the  Taylor  series  expansion,  (1 1)  becomes, 

'  Y  (b  =  0;  c  =  0)  =  -  A  ^cos  (kx  -  cut)  +  £  Ak  cos  (2kx  -  2cot)j  +  0(A3k2),  (13) 

Compared  with  (5),  the  set-down  effect  in  the  water-surface  profile,  1),  (the  second  term  on  the 
right-hand  side  of  (5))  does  not  appear  in  the  shoreline  profile  along  the  beach  surface  in  (13); 
instead,  the  shoreline  profile  contains  its  second-harmonic  component  Hence,  the  shoreline 
profile  has  features  of  the  peaked  runup  and  flattened  rundown.  Figures  2  a  and  b  demonstrate  the 
difference  between  the  water-surface  profile  at  the  initial  shoreline  (5)  and  the  shoreline  profile 
along  the  beach  surface  (13). 

Taking  ti  (at  the  shoreline)  =  -  Y  sin  (3  and  solving  (5)  for  Y  also  yield  (13).  Hence,  (13)  is 
consistent  with  (5).  The  second-harmonic  component  in  the  shoreline  profile  is  due  to  the 
interference  of  the  exponential  offshore  profile  with  the  linearly  sloping  beach  surface.  This  effect 
is  demonstrated  in  Fig.  3.  At  Y  =  0  (the  shoreline  location  for  the  quiescent  state),  (5)  predicts  that 
the  vertical  water-surface  elevation  is  sinusoidal  in  x  and  L  However,  the  value  of  T]  cannot  be 
negative  at  Y  =  0,  hence  during  the  rundown  phase,  the  shoreline  is  determined  at  the  intersection 
of  the  exponential  decay  of  the  profile  and  the  linear  beach  surface.  During  the  runup  phase,  the 
shoreline  is  determined  by  the  exponential  extrapolation  from  Y  =  0.  The  exponential 
extrapolation/interpolation  is  the  one  that  causes  the  features  of  the  second  harmonic  in  the 
shoreline  profile.  The  set-down  effect  in  (5)  occurs  in  order  to  satisfy  the  conservation  of  mass. 
This  result  may  be  important  when  field  data  of  shoreline  variations  are  analyzed  for  edge  waves. 

The  work  for  this  paper  was  supported  by  the  Office  of  Naval  Research  (N00014-87-K-0815). 
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Figure  1.  --  Schematic  drawing  of  the  fluid  domain. 

Figure  2.  --  Edge-wave  profiles:  a)  water-surface  elevation  at  y  =  0  based  on  (5);  b)  shoreline 
profile  based  on  (13).  Ak  =  0.5  for  exaggeration  of  the  nonlinear  effects. 

Figure  3.  --  Maximum  and  minimum  offshore  profiles  of  the  edge  wave.  Ak  =  0.5  for 
exaggeration  of  the  nonlinear  effects. 
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